Spontaneous emission by rotating objects: A scattering approach 
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We study the quantum electrodynamics (QED) vacuum in the presence of a body rotating along 
its axis of symmetry and show that the object spontaneously emits energy if it is lossy. The radiated 
power is expressed as a general trace formula solely in terms of the scattering matrix, making an 
explicit connection to the conjecture of Zel'dovich [JETP Lett. 14, 180 (1971)] on rotating objects. 
We further show that a rotating body drags along nearby objects while making them spin parallel 
to its own rotation axis. 
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Quantum zero-point fluctuations of the electromag- 
netic field in vacuum lead to macroscopic manifestations 
such as the Casimir attraction between neutral conduc- 
tors When objects are set in motion, they may pull 
out real photons from the fluctuating QED vacuum. In 
fact, accelerating boundaries radiate energy, and thus 
experience friction, through the dynamical Casimir ef- 
fect (see Ref. for a recent review) . Even two paral- 
lel plates moving laterally at a constant speed experience 
a non-contact frictional force [H, @j- While a constant 
translational motion requires at least two bodies (other- 
wise, trivial due to Lorentz symmetry), a single spinning 
object can experience friction. In a recent work [6|, Man- 
javacas and Garcia studied such rotational friction by ex- 
pressing the polarization fluctuations of a small spinning 
particle via the fluctuation-dissipation theorem, and ob- 
tained a frictional force even at zero temperature. In fact, 
this problem is closely related to a classicalphenomenon 
known as superradiance due to Zel'dovich [7j. He argues 
that a rotating object amplifies certain incident waves, 
and further conjectures that, when quantum mechanics 
is considered, the object should spontaneously emit radi- 
ation only for these so-called superradiating modes. In- 
deed this is shown to be the case for a rotating (Kerr) 
black hole by Unruh Q. This phenomenon, however, is 
different in nature from Hawking radiation Q. One can 
also find similar effects for a superfluid (lo| . 

In this letter, we treat the vacuum fluctuations in the 
presence of a rotating object exactly, except for the as- 
sumption of small velocities to avoid complications of rel- 
ativity. By incor por ating scattering techniques into the 
Rytov formalism ll[ , we find a general trace formula for 
the spontaneous emission by an arbitrary spinning ob- 
ject, solely in terms of its scattering matrix. We repro- 
duce (and correct) the results in the literature and find 
an expression for the radiation by a rotating cylinder. Fi- 
nally, we study the interaction of a rotating body with a 
test object nearby and show that the rotating body drags 
along nearby objects while making them rotate parallel 
to its own rotation axis. 

Our starting point is the Rytov formalism [ill ] which 
relates fluctuations of the electromagnetic (EM) field 
to fluctuating sources within the material bodies, and 



in turn to the material's dispersive properties, via the 
fluctuation-dissipation theorem. The EM fields are gov- 
erned by the Maxwell equation 



V x V x 



e(w,x)I 



(1) 



where a linear, and nonmagnetic medium is assumed. 
Rytov postulates that the sources undergo fluctuations, 
which are related to the imaginary part of the local di- 
electric response e(w,x) by 

(K(w,x)®K*(w,y)> =or(w)Ime(w,x)«(x-y)I (2) 



4h(n,T(uj) + 1/2) where nr(w) = 
_1 is the Bose- Einstein occupation 



Here, ar(w)) = 
[exp(Sw/fc B T) - 1 
number, and T is the temperature. Equations ([1]) and 
([2]) define the fluctuations of the EM-fields in the pres- 
ence of a static object at temperature T. For bodies in 
uniform motion, they are applied in the rest frame of the 
object, and then transformed to describe the EM-field 
fluctuations in the appropriate laboratory frame. With 
all contributions of the field correlation functions in a sin- 
gle frame, one can then compute various physical quan- 
tities of interest, such as forces, or energy transfer from 
one object to another, or to the vacuum. For non-uniform 
motion, we assume that the same equations apply locally 
to the instantaneous rest frame of the body [12j. This 
assumption should be valid as long as the rate of acceler- 
ation is less than typical internal frequencies characteriz- 
ing the object, which are normally quite large. The mod- 
ified Maxwell equations are easier to derive from a La- 
grangian, £em = h e' E' 2 — | B' 2 , with the primed fields 
defined with respect to the co-moving reference frame. 
To incorporate the fluctuating sources, we must add 



AC = K' • E' 



(3) 



to the Lagrangian, where K' is defined in the object's 
frame. To make contact with the EM-field in the vacuum, 
we should recast all fields in the stationary (lab) frame. 

In this letter we assume that the velocity of the ob- 
ject is small such that, to the lowest order in v/c, the 
electromagnetic fields transform as 



E' 



xB, B' = B 



x E. 
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We discuss below the limits that this assumption places 
on the generality of our analysis. In the lab frame the 
EM equation is then given by 



V x V x — -I - — D(e' - 1) 



E = ^DK', (4) 



where D = I + J-v x V x and = 1 + —V x vx. 
Also e' is the response function defined in the object's 
rest frame, now written in terms of laboratory coordi- 
nates. Equation Q describes the EM-field in the lab 
frame in terms of sources defined in the object's frame, 
(K'(u',x')®K'*((j',y')) =aT(w')Inie(w',x ; )^(x'-y')I, 
where uj' and x' are understood as the frequency and po- 
sition in the moving frame, which should be transformed 
to those in the lab frame. 

Next we turn to the main point of this study, namely 
a solid of revolution spinning with angular frequency O 
along its axis of symmetry. We choose time and polar 
coordinates (t, r, </>, z) and {t 1 , r', 4>', z') in the lab and the 
object reference frames respectively, and take the rota- 
tion along the z axis. The two coordinate systems are 
related by 



t' = t. 



r, (j)' 



fit, 



(5) 



Consider a fluctuation of the source characterized by 
frequency u/ and azimuthal index m', K^, m , (i',x') = 

e- iw ' t '+ im '*'f u ,',m'(r / ,z / )- Note that f depends only on 
the coordinates r' and z' . As a matter of notation, we de- 
fine K(t, x) = K'(i', x'), i.e. we drop the prime when the 
function is expressed in the lab-frame coordinates. The 
coordinate transformation in Eq. ([5]) then implies that 
the frequency as seen by the rotating object is shifted by 

K w . m (i,x) = K LJ _f2 Tn m (t ,x ). 

This equation, in turn, recasts the source fluctuations 
into the lab coordinates: 

(K m (cj,x) <g> K* m (uj,y)) = a T (u} - flm)x 

xlmetco- flm, r, z) S ^ x ~ ry) ^ x ~ Zy) I. (6) 

2irr 

Here, we used the constraint that the object is rotation- 
ally symmetric in the c/> direction. Equation (0| together 
with the above equation describes the electromagnetic 
field fluctuations in the presence of a rotating body. Note 
that Eqs. ^ and © include both positive and negative 
frequencies; since Ot(w) and Ime(u;) are odd functions, 
rest frame fluctuations are identical for opposite signs. 

The electromagnetic fields outside the object receive 
contributions both from the fluctuating sources within 
the object and from fluctuations (both zero-point and at 
finite temperature, thermal) in the vacuum outside the 
object. The source fluctuations in the vacuum are given 
by Eq. © as described below. 



First we consider the fluctuations inside the object and 
find the corresponding field correlation functions. Equa- 
tion (Ql together with free Maxwell equation in the vac- 
uum (e = 1, K' = 0) gives the electric field via the 
Green's function 



Ein-fluc(w, x) 



dzG(w,x,z) -BK(ij,z), (7) 



where we have rewritten K' in terms of lab-frame coor- 
dinates. The subscript on E indicates that the electric 
field is due to the inside fluctuations (but is possibly com- 
puted outside the object). The required Green's function 
(with one point inside and the other outside the object) 
can be formally expanded as 



C a 



G( W ,x,z)=£HrE° u >,x 



,(w,z) 



(8) 



The index a m denotes a set of quantum numbers includ- 
ing m, the eigenvalue of the angular momentum along 
the z direction (in units of h) . For example, for spherical 
waves a m = (P, I, m), where P is the polarization and I is 
the eigenvalue of the total an gula r momentum. The pref- 
actors Ca m as defined in Ref. [13( are introduced for later 
convenience. Here, E out ( m ) is the usual outgoing (incom- 
ing) wave, while F is a solution to the homogeneous EM 
equation inside the object, i.e. Eq. Q with the RHS 
set to zero. The important constraints are the continu- 
ity equations these functions satisfy on the surface of the 
object 



(F a J|| = (E' Q n m +S Qm Er„ 



where S am is the scattering matrix. A similar equation 
holds for the curl acting on the fields. It turns out that 
the knowledge of the surface values of F, which in turn 
can be expressed in terms of the scattering matrix S am , is 
sufficient for computing the correlation functions. Equa- 
tions (|7|) and (JSJ along with the source fluctuations of 
Eq. ([B]), and a few integrations by parts yield 

iC uJ 2 

(E(w,x) ® E*(w,y))i n _fi uc = ^2a T (uj - fim) — |^ — 



i-|s a J 2 )C( w . x )®CKy)' 



(9) 



Next we turn to fluctuations caused by outside (vac- 
uum) sources. Equation ([5]) may appear to suggest that 
these fluctuations vanish because Ime = in the vac- 
uum, whereas vacuum fluctuations exist even in the ab- 
sence of any objects. The key is that an integral over in- 
finite space leads to 1/ Im e and thus the limit Im e — > 



should be taken with care |14| . A careful analysis similar 



to that for the inside contribution yields the correlation 



3 



function due to the outside sources as 



(E(w, x) ® E*(w, y))out-fluc = dTn (w) 



4c 2 

x (E^( W ,x) + ^ m E™ t ( W ,x)) 
®(E 1 a -(a;,y)+S: m E- t *(a;,y)). (10) 

Note that the function aT is defined at the environment 
temperature To and depends only on w, unlike the inside 
fluctuations which depend on the shifted frequency, lo — 
Sim. Equation (|10l) has a rather intuitive form: each 
term in the dyadic expansion is a linear combination of 
the incoming wave plus the scattered outgoing wave. 

Summing Eqs. © and (fTU)) . the total correlation func- 
tion is given by 



(E® E* 



(E ® E* 



in- flue 



(E ® E*) out _ nuc . 



The above equation completely characterizes the vacuum 
fluctuations in the presence of a rotating body possibly 
at a different temperature from the vacuum. Interest- 
ingly, even at zero temperature, the rotating body spon- 
taneously emits energy, as can be computed by averaging 
over the Poynting vector, (E x B), and integrating over 
a surface enclosing the object. The sum of the in- and 
out-fluctuation contributions yields 



P 



2^ 



J^Mw-fimJ-nibCw)] (l-|S Qm ( W )| 2 ). (11) 

Note that the singularity of rir at u> = Sim is removed 
since 1— IS^ 2 vanishes there. Equation (fTTj) can be written 
in a basis-independent form, where the radiation takes 
the form 



P 



2^ 



hjj x 



Tr 



(n T (uj - Sll z ) - titoM) (I - S f (u;)§H) 



(12) 



Here l z is the z-componcnt of the angular momentum 
operator (in units of h) and S is the (basis-independent) 
scattering matrix. Note that this equation reduces to the 
heat radiation from a static object in the limit of zero 
angular velocity 3 HI]. We are specifically interested 
in zero temperature. In this limit, the number of radiated 
photons in the partial wave a m is 



duj 



= G(Slm - lo) {\S am (uj)\ 2 -1), 



(13) 



where is the Heaviside function. Hence, the radiation 
only comes from the frequency window < ui < Sim for 
the m th partial wave. In fact, Zel'dovich proposed classi- 
cal superradiancc for the same frequency regime (tJ. He 
considered a rotating cylinder and argued that exactly 



for waves in the above frequency range, the amplitude 
of the scattered wave (in absolute value )is larger than 
one, i.e. \S am (u)\ > 1; see Refs. [3, 17 for further 
discussion. He further conjectured that taking quantum 
mechanics into account would lead to spontaneous emis- 
sion. Equations (|1HI13[) are, to our knowledge, the first 
expressions for spontaneous emission by a rotating ob- 
ject given explicitly in terms of the scattering matrix, 
which show that, at zero temperature, the radiation is 
generated exactly in the superradiating channels. 

Note that, having assumed non-relativistic velocities, 
we must keep only the leading contribution in SIR/c to 
radiation from each partial wave. In fact, higher par- 
tial waves typically contribute in higher powers of this 
quantity. Therefore, to compute the total radiation at 
zero temperature, we must keep only the lowest partial 
wave, while at finite temperature higher partial waves 
can make a comparable, or even larger, contribution due 
to the Boltzmann weight. 

Using the general expression for radiation, we now dis- 
cuss some simple special cases, namely a sphere and a 
cylinder. To find the 5-matrix of a rotating object we 
have to solve a complicated equation — Eq. (j4]) with the 
RHS set to zero, but the task is made easier by mak- 
ing a further assumption that the object's radius is small 
enough so that \y/e\SlR/c -C 1, which allows us to neglect 
the explicit dependence on the velocity (D w I). Also, 
for small ft we need only consider the low-frequency re- 
sponse. 

For a sphere the T-matrix is related to its (electric) 
polarizability a by: T® n ® lm {uj) = i^-a(oj — Sim): the 
argument of a is u> — Sim because that of e is shifted. 
The value of these functions for negative arguments can 
be obtained through their analytic properties. The scat- 
tering matrix is related to the T-matrix via 5 = 1 + 2T, 
and the energy radiation in the lowest partial wave can be 
computed by Eq. (jlip . The result is indeed in agreement 
with Ref. @. 

The T-matrix for the cylinder is rather complicated. 
Further, one must take into account all polarizations. 
Here, we quote the final result for a slowly rotating cylin- 
der of radius R and length L, and at zero temperature 
everywhere, 



P 



4HLR 2 
3irc 3 



duj uj 



Im 



e(w - SI) - 1 



e(u) - ft) + 1 



(14) 



This equation, valid for arbitrary e, reproduces the result 
in Ref. [l8| in the limit of small conductivity. (Zel'dovich, 
et.al. miss a factor of one-half in the last step to their 
Eq. (3.7).) 

The radiation from a rotating object exerts pressure 
on nearby objects. Consider a spherical object with an- 
gular velocity ft along the z direction and a second small 
spherical body — a test object — at rest, placed at a sep- 
aration d on the x axis. One finds that, in addition to the 



4 



force along the x axis, there is a tangential force in the 
y direction. Furthermore, a torque is exerted on the test 
object. To compute this effect, we break up the EM-field 
correlation function into radiation (due to propagating 
photons) and non-radiation (due to zero-point fluctua- 
tions) parts: 



(E <g> E* 



(E ® E*) ra d + (E®E* 



1 non— rad • 



The non-radiation part gives the Casimir force which 
makes no contribution to the tangential force or the 
torque. At zero temperature, the radiation comes only 
through the superradiating modes with uj < f2m, 



(E( U ,x)«E*( W ,y)) rad = fi^e(fim- [ d) 

OCm 

x (1 - \S a J 2 ) E° u >,x) ® E°f ( W ,y). 



(15) 



We assume that d is large compared to the length scales 
of both the rotating body and the test object, and thus 
compute the first reflection of the radiation off of the 
test object [19j . First we transform the outgoing waves 
to regular waves about the origin of the second object 
via translation matrices [13( . These waves scatter on the 
second object giving 



M 



he 2 
8ird 2 



1 



duj—(\SnE\ 2 -l)(l-\&iiE\ z ), (16) 



for the torque, and 



32-Kd 



/ cM|Sii £ | 2 -l)(l- Re6ii£), (17) 
Jo 



for the shear force on the test object, where She an d 
Si ib are the scattering matrices of the rotating and the 
test object, respectively, in the lowest partial wave (Z = 
m = 1, P = E). This partial wave gives the leading order 
in 1/d for small f2. Note that the torque falls off as 1/d 2 
with separation while the force goes as 1/d. Furthermore 
the signs are positive, i.e. a rotating body drags along 
objects nearby while making them rotate parallel to its 
own rotation axis. 

To get an estimate for the magnitude of radiation ef- 
fects, we consider a rapidly spinning nanotube of radius 
R and length L, and assume that CIR/c is small. We 
then find that the rotation slows down by an order of 
magnitude over a time scale of t ~ {I /ft) (c 3 / LR 2 il 3 ). 
The moment of inertia of a nanotube can be as small as 
1(T 33 in SI units [U (compare with h w 1CT 34 ). So even 
at small velocities, r can be of the order of a few hours. 

We have derived a universal formula for the sponta- 
neous emission by a rotating object in terms of the scat- 
tering matrix that makes an explicit connection to the 
physical principles behind supcrradiancc. Furthermore, 
it allows one to circumvent assumptions about small size 



or small conductivity used in the literature. We also be- 
lieve the current formalism naturally generalizes to rel- 
ativistic motion. Finally, generalization of the technical 
and conceptual aspects of this work to the interaction of 
multiple (moving) objects would be worthwhile. 
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